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Abstract 

By using the differential operator technique and the effective field theory scheme we 
study the tricritical behavior of Heisenberg classical model of spin-1/2 in a random 
field. The phase diagram in the T-h plane on a square and simple cubic lattice for a 
cluster with two spins is obtained when the random field is bimodal distributed. 
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1. Introduction 

The random field Ising model (RFIM) has been investigated extensively both theoretically 
and experimentally in last yearsS ~§ Some questions such as the lower critical dimension 
and the existence of a static phase transition have been already solved from the theoretical 
point of view, but others questions as the existence of the tricritical point are still open. 
However, it is now understood that RFIM describes the essential physics of a rich class of 
experimentally accessible disordered systems. These include structural phase transitions in 
random alloys, commensurate charge- density- wave systems with impurity pinning, binary 
fluid mixtures in random porous media, and the melting of intercalates im layered compounds 
such as TiS2. The behavior of the RFIM is governed by the competing tendencies of the 
spins to align ferromagnetically under impetus of the exchange or to follow local fields and 
so be uncorrelated. 

The RFIM has been studied with different methods. The renormalization-group the- 
ory (RG) of phase transition and critical phenomena represent one of most powerful tool 
of theoretical physics in last decades. Until then (early 1970s) all calculations of critical 
exponents were made either by exactly or by series expansions. Therefore, employing the 
renormalization-group arguments (which can be applied for weak fields), the proof of the 
equivalence between an Ising ferromagnet in a random field and a dilute antiferromagnet in 
a uniform field has been demonstrated. These question have been very explored (see |] and 
references therein). 

In the RG scheme, such as, the effective field RG (EFRG)ii and the Mean-field RG 
(MFRG)i'i methods have been successfully employed in spin systems, where the results 
obtained are in accordance with other more effective approach (series expansion, Monte Carlo 
simulation among others). Those approaches (MFRG and EFRG) are based on comparison 
of two clusters of different sizes, each of them simulating infinite systems. On the other 
hand, the mean field approximation (MFA) has been extensively applied in the study of 
RFIM. The MFA provides very interesting results about RFIM. For example, a bimodal 
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distribution for random field gives rise to a tricritical point and consequently first order 
phase transitions, at very low temperatures for certain values of field strength, but does not 
occur in the case of a gaussian distribution.0~0 In addition, it is possible to establish a 
complete mapping between the parameters of the Ising ferromagnet in a random field and 
the dilute Ising antiferromagnet in a uniform field. 

By employing the same strategy of EFRG, the effective field theory (EFT) has been used 
to study cooperative phenomena and phase transitions in various systems and giving useful 
qualitative and quantitative insights for the critical behavior of a wide variety of classical 
and quantum lattice spin systems.lli~0 The approach of EFT uses, as a starting point, the 
rigorous Callen-Suzuki spin identities (See Refs. |18|and|l^ and the effects of the surrounding 
spins on each cluster are taken into account by using a convenient differential operator 
technique introduced by Honmura and Kaneyoshi.0 In this procedure, all relevant self-spin 
correlations are taken exactly into account in the EFT equation of state.ill Therefore, the 
EFT approach is superior to the standard MFA. 

In this work, we study the behavior of phase diagram {in T — h plane) for Heisenberg 
classical model of spin-1/2 on a square and simple cubic lattice by employing a bimodal 
probability distribution for random fields in a cluster with two spins by using the result 
recently obtained on diluted systems.l^ The outline of the remainder of this paper is as 
following: in section 1, the formalism and calculations are developed and the results and 
conclusions are presented in section 2. 

2. Formalism and Calculations 

The system to be studied is the n- vectorial model in a random field described by the Hamil- 
tonian 

-/3H = i^^Si-Sj+^hi-Si, (1) 

(hi) » 
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where the summation is performed over all pairs of the nearest-neighboring sites The 
quantities Si are isotropically interacting n -dimensional classical spins of magnitude y/n 
localized on site i and the Cartesian components of Si obeys the normalization condition,ii 



n 



n\2 



n. 



K [= J/ksT, ks is the Boltzmann constant and T the temperature] is the exchange in- 
teraction between the spins, hi [h = fisH/ksT , where fi is the Bohr magneton and H is 
the random magnetic field] is the reduced random magnetic field at site i with probability 
distribution 

P{hi) = ^[S{h, + h)+6ih,-h)]. (2) 

Hamiltonian (|l]) reduces to the S = ^ Ising, planar (XY), Heisenberg and spherical 
models for n = 1, 2, 3 and oo, respectively. 

In this work, we follow the EFT procedure (see Ref. p!6| ) to study the critical properties 
of the Hamiltonian described by Eq.(|l]) by employing the axial approximation.^ Since the 
Hamiltonian for a cluster with two spins can be written as 

H = KSi-S2 + ai- Sl + a2- Sl, (3) 

where ai = hi + K^^jS], (/ = 1, 2) and z is the lattice coordination number. 

The average magnetization per-spin m = ^{{Sl + SD) for cluster with two spins is given by 

(See Refs. ^ and |^) 

^ = ( n + S'A) n + Slf3y))gniX, Y)\x=huY=h,. (4) 

= cosh{KD^), (3^ = smh{KD^), {u = x,y), X{Y) = x{y) + hi{h2). {= is 
differential operator which satisfy the mathematical relation 

sinh{aD^ + hDy)gn{X,Y)\x=h^,Y=h2 = 9n{a + hi,b + /ig), 
where gn{X, Y) is given by 



a(XY) = ^i^M^ + Y) 

' ^ cosh(X + F)+exp(-2ir)M„(ir)cosh(X-F)' 

and 



ln/2{nK) - In/2-iinK) 



In/2{nK) + In/2-l{nK) 

Here is a modified Bessel function of the first kind. Eq. (^) is exact and will be 

applied here as a basis of our formalism, since it yields the cluster magnetization m and the 
corresponding multi-spin correlation functions associated with various sites for the cluster 
under consideration. Here we apply the EFT approximation in both sides of Eq. (P, i.e., 
the thermal and random average (denoted by (■ • ■)c)^ along with the decoupling procedure 
which ignores all high-order spin correlations, namely {SjS^ ■ ■ ■ S^c ~ {S})c{Sj)c ■ ■ ■ {Sn)c , 
with i ^ j ^ ■ ■ ■ ^ n. Based on this approximation and replacing each boundary configu- 
rational spin average by the symmetry breaking mean-field parameters 6j for all i, one set 
up the equation of state for rh = (m)^. By using the properties of differential operator and 
assuming translational invariance, we obtain for the square (sq) and simple cubic lattice 
(sc), respectivelly, 
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rh = '^Ak,sq{K,n,h)Th'', (6) 

k=l 
10 

m = Y,Ak,sc{K,n,h)Th\ (7) 

k=l 

The coefficients Ak,(sq,sc){K,n, h) [k = even] are zero to satisfy the time reversal symmetry 
of the Ising model, as well the properties of differential operator technique. 
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3. Results and Conclusions 



In this section we apply the conditions for determining of the second-order transition hne 
and the tricritical point (TCP) in the n— vector model in a random field. We will focus on 
the case n = 3, which corresponds to the Heisenberg classical model of spin-1/2. 
In the vicinity of the second-order phase transition, the Eqs. (J^ - |^ is given by 
^2 _ h) - 1 

M,{sq,sc){K,n,h) 

Since the magnetization (m) goes to zero continuously, to obtain the second order phase 
transition line, in the T — h plane, we need solve the equations 

Ai^(sq,sc){K,n,h) = 1, A-i^^sq,sc){K,n,h) < 0. (9) 

On the other hand, the r.h.s. of Eq. (|^) must be positive; otherwise, the transition is 
interpreted to be of first-order. Unfortunately the first order transition line cannot to be 
analyzed in the frame of the differential operator technique, therefore we have confined our 
calculations only to the second order transitions, including the TCP. We obtain the TCP 
by solving the equations 

Ai^^sq,sc){K,n,h) = 1, A3^^sq,sc){K,n,h) = 0. (10) 

The numerical solution of Eq. (H) provide the second-order phase transition line which 
is shown in Figure 1. When the random field is zero, we obtain the well-known critical 
temperature of Heisenberg classical model of spin-1/2, K^^^ = 3.005 and K~}^ = 5.031. 
In Figure 1, the solid line indicate the second-order phase transitions; the black diamonds 
denote the position of the TCP at which phase transition changes from second to first 
order. The numerical solution of Equations (|1^) provide the values for the tricritical points 
{hTCP,KTCp) = (1.800,0.958) and {hTCP^K^hp) = (2.785,2.389) for the square and simple 
cubic lattice, respectivelly. We might conclude saying that EFT approximation is able to 
predict the presence of tricritical behavior for the Heisenberg classical model of spin-1/2 on a 
square and simple cubic lattice in a random magnetic field which takes on the random values 
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±h with equal probabilities. As can be seen, our results show the existence of the TCP for 
square lattice {z = 4) in accordance with mean- fieldli and Betthe- Peierlsil approximat ions . 
This occur because to obtain the first-order phase transition we must have at least the rh^ 
term in the expansion of Eqs.(^- Therefore, in the presente approach we cannot observe 
behavior tricritical for lattice with coordination numbers z < 4. 
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Figure Captions: 



Figure 1. Phase diagram for Heisenberg classical model in presence of the random mag- 
netic field, in the T — h plane, on a square (sq) and simple cubic (sc) lattice. 
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